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Рассматривается математическое моделирование и идентификация параметров динамической связи между 
инструментом и обрабатываемой заготовкой, формируемой в процессе резания. Приводятся данные о свой- 
ствах динамической связи при малых вариациях координат в окрестности точки равновесия. В статье рас- 
сматривается исключительно позиционная связь. 
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Введение. Математическое моделирование динамической связи, формируемой процессом реза- 
ния и объединяющей динамические подсистемы станка со стороны режущего инструмента и об- 
рабатываемой заготовки, являлось предметом многочисленных исследований [1-4]. Это связано с 
тем, что формируемая процессом резания динамическая связь во многом определяет устойчи- 
вость траекторий стационарных формообразующих движений инструмента относительно заготов- 
ки, и она влияет на преобразование траекторий исполнительных элементов станка в траектории 
формообразующих движений. Для изучения динамики процесса резания используется следующее 
обобщенное уравнение: 


2 
Х х . 
пе -+о(Х) = РУ», б,иь)+ У (6), (1) 
где Р(Х,Г,, 5,1) = Е (ХИ, 5), ВХ, ь 1)... ВЫХ ,Гь,5ь,»)} — вектор-функции ди- 


намической характеристики процесса резания, раскрывающие зависимость сил резания от 
упругих деформационных смещений инструмента и заготовки, а также от технологического 


режима: величины подачи на оборот 5, глубины резания {»и скорости резания Г» ; 
ХЕХ А, а — вектор упругих деформационных смещений вершины режущего ин- 
струмента (первые три координаты) и заготовки в точке контакта с ней режу-щего инстру- 
мента (последние три координаты); /“(#) ={/“(#), №/(0,..., /@)' - изменяющиеся во 
времени составляющие сил резания, не объяснимые в координатах упругих деформацион- 
ных смещений, которые интерпретируются как шум; т(Х) = [2 2] #(Х) = [и (Х)| 
с(Х) = [с (Х)| 5, К =1,2.....6 — соответственно функциональные матрицы инерционных и 


диссипативных коэффициентов, а также функциональная матрица формирования упругой 

составляющей сил в зависимости от вектора деформационных смещений, которые, в свою 
Т 

очередь, зависят от внешних сил. Во всех случаях здесь и ниже символ {...}? - есть опера- 


ция транспонирования. 
Вопрос о математическом моделировании и идентификации функциональных матриц 


т(Х) = [А | (Х) = [и (Х)| с(Х) = [с (Х)| рассмотрен нами ранее [5, 6]. В настоя- 
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щей статье рассматривается проблема математического моделирования и идентификации вектор- 
функций Е(Х,Г,,5,,1»). При этом используются алгоритмы скользящей линеаризации, рас- 


сматривающие линеаризованные их оценки в окрестности точки равновесия д. задаваемой 
технологическими режимами. Исследования построены таким образом, что в рассматриваемом 
частотном диапазоне параметры жёсткости подсистемы заготовки на порядок превышают жёст- 


кость подсистемы режущего инструмента. Поэтому уравнение (1) в вариациях относительно Х“ 
(Х@=Х" +х(@) ) можно представить в линеаризованном виде 

ах 
а!” 
где х=(х,,х,,х.,)’ — вектор упругих деформационных смещений вершины режущего инструмен- 





* * [о * 
т АА, (Х „Ир бр)+сь(Х И, ,5 5 1ь)х = ГВ), (2) 


та В вариациях относительно ТОЧКИ равновесия, причём Х* = с0п% р 


т(Х”)= | (х')| 5, К =1,2,3 -матрица масс подсистемы инструмента; 


ВХ") = (ХИ, бь >) 5. =1.2,3, © (Х") = в", бьь) |8, = 12,3 = 
суммарные матрицы коэффициентов демпфирования и жёсткости, состоящие из параметров 
подсистемы инструмента и параметров динамической жёсткости и демпфирования процесса 
резания; /(г) ={/ (0) Х, (#), Г. (Е) - вектор внешних измеримых дельтообразных силовых 
возмущений. 
В (2) суммарные матрицы: 
(ХР, вь) = ИСХ") + (ХИ, 5), 
мо а (3) 
где #“® = [и = [бФ, /ах, 5,к=1,2,3, с = [|= [0Ф,/дх,|] $,к=1,2,3 -— соответственно 
матрицы скоростных коэффициентов и жёсткости процесса резания, рассматриваемые в ва- 


риациях относительно точки равновесия. 
Следовательно, 
ф.(х) = РЁ. (И, 51», К +х)- Е, (Ть, 51», Х"), 5=1,2,3 и ф, (0) =0. 

В статье рассматриваются свойства динамической связи, формируемой процессом реза- 
ния, в линеаризованном представлении для различных точек равновесия, задаваемых технологи- 
ческими режимами, - в этом смысл скользящей линеаризации, справедливость которой определя- 
ется тем, что при динамических исследованиях смещение точек равновесия есть медленная 
функция времени. В рассматриваемом частотном диапазоне, ограниченном сверху частотой по- 
рядка 5,0 кГц, инерционные реакции со стороны процесса резания не наблюдаются. Уточним по- 


нятие точки равновесия Х” динамической системы резания. Это точка, в которой скорости и ус- 
корения упругих деформационных смещений вершины инструмента равна нулю. Этой точке соот- 
ветствует движение обрабатываемой заготовки с постоянными скоростями резания и подачи. Ста- 


тья посвящена исключительно изучению с“” =[® = [9Ф, дк | 5,к =1,2,3. Идентификация и 


свойства матрицы скоростных коэффициентов будут рассмотрены в следующих наших публикаци- 
ях. 
Идентификация и свойства матрицы динамической жёсткости. Идентификация элементов 


матрицы с“”(Х”,Г,,5 1») осуществляется на основе анализа реакций системы в процессе ре- 
зания на внешние силовые возмущения, в том числе, дельтообразные силовые возмущения. В по- 


следнем случае рассматривается три совокупности сил: /({) = {1 (#),0,01', Х(®) = {0, 1. (#),0' 
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и Г(#) = 0,0, Г, (#)}', которым соответствуют три совокупности вибрационных реакций инстру- 
мента, рассматриваемых в виде колебательных ускорений, измеренных с помощью пьезоакселе- 


рометров. Это следующие совокупности: х‘? (9) = о () 15, К=1,2,3 , которые близки ко вторым 


производным от весовых функций. Кроме этого анализируются их частотные образы в виде спек- 
тральных матриц колебательных реакций и также сил и колебаний. Идентификация выполняется 
в два этапа. На первом этапе определяются матрицы т, Й, с в системе без резания. Методика и 
математическое моделирование этих матриц приведена в наших публикациях [5,6]. На втором 
этапе оцениваются матрицы с; и непосредственно матрицы динамической жёсткости процесса 


резания с; —с = с” . Таким образом, матрицы динамической жёсткости оцениваются на основе 


анализа весовых функций или их Фурье-образов в низкочастотной области, лежащей существенно 
ниже собственных частот подсистемы инструмента. Это частотный диапазон, в нашем случае ог- 
раниченный сверху частотами порядка (200-300) Гц. Полученные оценки сравнивались с отноше- 
ниями соответствующих составляющих сил резания к смещениям инструмента, обусловленным в 
установившемся состоянии соответствующими вариациями технологических режимов (величины 
подачи наоборот и глубины резания с учетом упругих деформаций). 

При оценивании матрицы динамической жёсткости перечислим основные факторы, 
влияющие на формирование сил резания. Они характеризуют следующие свойства, которые не- 
обходимо учитывать при построении динамических моделей процесса резания. 

1. На силы резания, прежде всего, оказывает влияние площадь срезаемого слоя. При рас- 
смотрении деформационных смещений и сил будем использовать системы координат, показанные 
на рис.1, а, 6. Тогда малые деформационные смещения инструмента по оси Х, не вызывают из- 
менения площади срезаемого слоя. В процессе точения обычно выполняется условие 


15/9 р)(10,0 — 20,0) . В этом случае чувствительность изменения сил к вариациям смещений ин- 


струмента относительно заготовки в плоскости (Х’, ‚ Х, ) принципиально зависит от направления 
деформационного смещения. Она максимальна по направлению нормали к режущей кромке 
инструмента и минимальна по направлению её ориентации (на рис.1 — это оси (Х®,Х®), 
повернутые на угол ф против часовой стрелки). На рис.1, в, гданы диаграммы приращения сил в 


зависимости от постоянного значения смещений инструмента в установившемся состоянии, рав- 
ных 0,1 мм, имеющих различные направления. Смещения инструмента обеспечиваются варьиро- 
ванием технологических режимов с учётом деформационных смещений инструмента. Во всех слу- 


чах ф=60°”. Приведены два примера. На рис.1, в даны диаграммы для соотношения 


1/5 р =3,3 , а на рис.1, г- для 1» /5› = 20,0. Там же точечными диаграммами приведены ап- 
проксимации изменения сил в виде эллипсов. Как видно, во втором случае, наиболее характер- 
ном для процесса точения, деформационные смещения в направлении оси Х“ практически не 


вызывают изменения сил. 

2. Математические модели подсистемы режущего инструмента построены таким образом, 
что параметры модели идентифицированы в частотном диапазоне, определяемом первыми фор- 
мами колебаний. 
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Рис.1. Схема оценивания вариаций сил и деформационных смещений в плоскости (ХЛ, ,Х’, } 


Однако при больших вылетах инструмента из зажимного приспособления приходится счи- 
таться с изгибными деформациями тела инструмента, сопутствующими линейным деформацион- 
ным смещениям. В этом случае деформационные смещения дополнительно вызывают изменения 
углов резания. Здесь прежде всего, необходимо учитывать изменения переднего угла инструмен- 
та, заметно влияющего на модуль и ориентацию силы. Причём зависимость сил от изгибных де- 
формационных смещений является нелинейной. Более того, при значительных изгибных дефор- 
мациях увеличению деформационных смещений может соответствовать не уменьшение сил, а их 
увеличение. Однако в вариациях относительно точки равновесия все характеристики вариаций 
сил, вызванные малыми деформационными смещениями можно линеаризовать. 

Таким образом, в зависимости от технологического режима и формы колебаний меняется 
структура матрицы динамической жёсткости процесса резания, представляющая собой коэффи- 
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циенты влияния деформационных смещений на составляющие сил различного направления. В 
общем случае матрица динамической жёсткости процесса резания представляется в виде 


со = [| [ф.о зК=123, (4) 
где столбцы матрицы (4) изменяются в зависимости от технологических режимов и форм дефор- 
мационных смещений. Здесь можно рассматривать три варианта. Первый относится к случаю, ко- 
гда выполняется условие {,/5,)(10,0—20,0) и изгибные деформационные смещения инстру- 
мента отсутствуют. Тогда в матрице (4) значимыми являются только коэффициенты первого 
столбца. Остальные равны нулю. Второй относится к случаю, когда 1, /5,((10,0—20,0) и ли- 
нейные деформации существуют вместе с изгибными. Тогда в (4) значимыми являются первый и 


третий столбцы. Наконец, наиболее общий случай, когда необходимо учитывать все столбцы мат- 
рицы (4). Он относится к случаю, когда дополнительно необходимо учитывать силовые реакции 


при деформационных смещениях в направлении Х:”. Самостоятельное значение имеют иссле- 


дования изменения матрицы (4) при варьировании скорости резания. 

Вначале рассмотрим матрицы динамической жёсткости для первого случая. При оценива- 
нии матриц динамической жёсткости рассматривается низкочастотная часть спектральных харак- 
теристик. Можно использовать также легкоизмеримые деформационные смещения вершины ин- 
струмента, оцениваемые с помощью тензодинамометра, показания которого тарируются в дефор- 
мационных смещениях вершины инструмента относительно станины станка. В этом случае де- 
формационные смещения необходимы для коррекции действительных значений величины подачи 
на оборот с учетом деформационных смещений. При идентификации матриц динамической жёст- 


кости  измеримыми являются силы  /“(0)={/(0),0,0', Г°(0) =, Г, (0),0" и 
1) (0) = {0,0, /,(0)}", которым соответствуют измеримые три совокупности деформационных 
смещений: 
(0 (0—9. (0 (0 т 
х (0 =ф: (0),х, (0), х5 (0)}, 
(2) (0) — 1.2) (2) (2) т 
х (0) =4х (0),х; (0), ху (0)}, 
(3) (()) — 5-8) (3) (3) т 
х (0) ={хг (0),х> (0), ху (0)}. 
Они рассматриваются в вариациях относительно точки равновесия, задаваемой технологическим 
режимом. Известными являются матрицы динамической жёсткости подсистемы инструмента с. 
При проведении экспериментальных исследований удобно анализировать не абсолютные значе- 
ния деформационных смещений, а их отношения. Это связано со сложностями обеспечения иден- 


ТИЧНОСТИ дельтообразных возмущений в ходе экспериментов. Кроме этого, при измерении коле- 
бательных ускорений их отношения равны отношению колебательных смещений. Тогда элементы 


матрицы с; можно вычислить из системы 


Ас=В, (5) 
ы т 
0 (5 а) = (5. 65) 
= : А= 9022) 0 — 922) 
где С= {Сих»бох»Сз;} А = ( 12 Сзз) ( 12 Сз1) , 
1,2 1,2 
(5. с) 0 — (65; 6) 


1 2 1 

Е 50 = м (0) (2,2) _ м (0) (1,2) _ г (0) | 

ПИЗЕН Е: = оз по" 
х› (0) х> (0) хз (0) 
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О Л (0) 
О) 


Подчеркнём, что в (5) нет параметров, в которых в явном виде выражаются силы или ко- 


т 
В = {(С›зСз — С›2Сзо ), (СозбСи эенен ро, > (Сб —С 46} го 


лебательные смещения. Возможность вычисления параметров {С .;,С\;,С,3х;} связана с тем, 


что элементы матрицы с известны. Кроме этого при составлении матрицы необходимо обеспечи- 
вать её хорошую обусловленность. Это зависит от параметров матрицы с и выбранных соотно- 
шений для относительных деформационных смещений. Например, если в качестве третьего урав- 


3 
(3,3) _ ж (0) 


нения взять соотношения для относительных деформационных смещений 0,5” = 3 (0) 
- х 


‚ то, как 


нетрудно показать, определитель от матрицы А всегда равен нулю. 

Параметры {с ;,С,›;›С, 3х} идентифицируются для конкретной точки равновесия, оп- 
ределяемой в установившемся состоянии технологическими режимами. В рассматриваемом случае 
при неизменной скорости резания — от ширины срезаемого слоя, то есть при неизменном угле ф 
— от величины припуска. Приведём пример изменения параметров {с ;,С,›;,С,з>} Для случая 
точения на станке 1Кб2 стали 20Х (рис. 2). Величина подачи на оборот во всех случаях равна 
$0) = 0,05мм/ об , скорость резания Г, =1.0м/с. Глубина резания варьировалась для полу- 


чения различных величин ширины срезаемого слоя. Геометрические параметры режущего инст- 
румента: нормальный передний угол (передний угол в нормальной секущей поверхности) 


0 № № ` 0 № 
7Тн и ‚ заднии угол в главной секущеи плоскости Ск =“. ; угол наклона режущеи кромки 


А =0°; главный угол в плане ф = 60°. Обрабатывается заготовка диаметром 4 =120мм ‚ за- 
креплённая в трехкулачковом патроне и поджатая вращающимся центром, после предваритель- 
ной обработки для устранения вариаций припуска. Приведём также данные по коэффициентам 
ориентации сил у, = Е [(Е)* +(Е,)* +(Е,)?]°°,1=12,3 (рис.3), соответствующие данным, 
приведённым на рис.2. 
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Рис.2. Пример изменения элементов матрицы жёсткости 
процесса резания при изменении ширины срезаемого слоя 
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Рис.3. Пример изменения коэффициентов ориентации сил при изменении ширины срезаемого слоя 


На приведённых иллюстрациях также даны точные значения математических ожиданий 
коэффициентов ориентации сил. Как видно, коэффициенты ориентации сил являются малоизме- 
няющимися функциями ширины срезаемого слоя, а зависимости элементов матрицы жёсткости от 
ширины срезаемого слоя в диапазоне изменения (1,0-5,0) мм можно линеаризовать. В последнем 
случае коэффициенты ориентации сил можно считать неизменными. Здесь вариации коэффици- 
ентов ориентации сил при изменении ширины срезаемого слоя не превышают 0,15 от номиналь- 
ного значения. 

В этом случае элементы матрицы жёсткости процесса резания при любом значении ширины 
срезаемого слоя можно задать величиной обобщённого коэффициента жёсткости и вектором ори- 
ентации сил. Для приведённых на рис.2 и 3 данных значение коэффициента жёсткости в диапа- 
зоне изменения ширины срезаемого слоя (1,0-4,0) мм можно представить в виде 


ео =рбх,1=1.2,3, (6) 

где р =190 - коэффициент, имеющий смысл давления стружки на переднюю поверхность инст- 
румента, кг/мм; Б - ширина срезаемого слоя, мм; у={7, 1..1.’ = {0,57;0,71;0,40}'. 

Подчеркнём, что коэффициент р’ зависит от скорости резания, геометрии инструмента и 


физико-механических свойств обрабатываемого материала. Заметим, что в процессе резания в 
связи с развитием износа инструмента указанные коэффициенты меняются. Это один из факторов 
эволюционного преобразования параметров динамической характеристики процесса резания. 


Кроме этого необходимо отметить, что угловые коэффициенты у ={7,,1.,х.}’ принципиально 


зависят от геометрии инструмента. В частности, по мере увеличения угла у„ =2° коэффициент 


Р К 
е не только уменьшается, но и может изменить свои знак. 


При оценивании матриц динамической жёсткости в общем случае необходимо восполь- 
зоваться прямыми измерениями изменения смещений и сил. Тогда для оценки суммарных значе- 
ний матрицы жёсткости можно воспользоваться системой 


ай (7) 
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оо бо Фобос оос 
9 0 0“ ^0 0.2 0 
0 0х’ 0 ох 0 0.57 
Хх 0 0 хо 0 0 хо 0 0 
Х= 0 хо 0 07 20 0:2 1 20 
0 0-я” 0 0-50 0 03? 
х" . 50 ох 0 о о 0 
0 0 ох 0 оо 
[0 0 хо 0 0 = 50 1. 
ных смещений; с= {Сила Сал Сола Соля Сол Ср Сзал>Сазт} = оцениваемые пара- 
метры суммарной матрицы динамической жёсткости; р ={/(0),0,0,0, }, (0),0,0,0, №. (0) = 
вектор внешних силовых возмущений, в частности, обусловленный изменениями технологи- 
ческих режимов. 

Проанализируем изменения матриц динамической жёсткости процесса резания при 
варьировании величины подачи на оборот (рис.4). Обработка осуществляется в условиях, соот- 
ветствующих рис.2, но глубина резания во всех случаях равна 1, =2,0 мм, что соответствует ши- 
рине срезаемого слоя Ь =2,5 мм. Варьировалась величина подачи на оборот, приводящая к уве- 


личению толщины срезаемого слоя а. На приведённых иллюстрациях сплошными линиями пока- 


заны коэффициенты Вий: = 1,2,3 , а пунктирными — коэффициенты в: = 12,3. 


где — матрица измеренных деформацион- 
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Рис.4. Пример изменения элементов матрицы жёсткости процесса резания при изменении величины 
подачи на оборот при различных значениях заднего угла инструмента 
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Обратим внимание на следующие особенности изменения элементов матрицы динамиче- 
ской жёсткости процесса резания. Во-первых, они принципиально зависят от заднего угла режу- 
щего инструмента. Во-вторых, по мере увеличения толщины срезаемого слоя, зависящей в основ- 


ном от величины подачи на оборот, имеет место существенное увеличение элементов матрицы 


Р Р Р 
Е и а" При этом элементы е) >) остаются практически неизменными. В третьих, при малых 


толщинах срезаемого слоя элементы матрицы в =1,2,3 практически равны нулю, и они не- 


пропорционально возрастают при увеличении толщины срезаемого слоя. 
В том случае, когда изгибная жёсткость инструмента становится существенной, приходит- 


ся считаться и с элементами матрицы, соответствующими второму столбцу. Для наглядности при- 


) Р) 


ведём зависимости изменения матриц с“” = [с | =1,2 для свободного резания на операции 


отрезки диска инструментом, имеющим специальную фаску в резцедержавке, способствующую 
значительным изгибным деформациям инструмента, вызывающим заметные изменения переднего 
угла инструмента (рис.5). 
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Рис.5. Изменение элементов матрицы динамической жёсткости процесса резания при отрезке диска 
при свободном точении в случае, когда изгибные деформации инструмента значительны. 


Точение стали 20Х при скорости резания 1,2 м/с и ширине срезаемого слоя Б = 3,0 мм 


На приведённой иллюстрации по мере увеличения сил резания, связанных с увеличением 


подачи в пределах 0,05-0,6 мм, передний угол ‘у уменьшался с 20° до 0. В результате принципи- 
ально изменяется структура матриц динамической жёсткости, зависящая от величины подачи на 
оборот. Во-первых, при малых подачах второй столбец матрицы с практически равен нулю, 
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что наблюдается в рассмотренных выше примерах. Однако по мере увеличения сил, связанных с 
увеличением подачи на оборот, значимым становится и второй столбец, что вызывает принципи- 
ально иные свойства динамической системы. Во-вторых, при этом меняется ориентация силы в 
пространстве, то есть меняются матрицы угловых коэффициентов её ориентации. Более того, при 
больших деформационных смещениях может измениться знак элементов матрицы динамической 
жёсткости. 

Выполненные исследования показали, что при больших изгибных колебаниях инструмента 


(Р) 


свойства матрицы с принципиально зависят от смещения точки равновесия в динамической 


системе. Как и в ранее рассмотренных примерах (см. рис.2, 4), так и практически при всех точках 
равновесия, задаваемых величиной подачи на оборот в примере рис.5, матрицы динамической 
жёсткости не обладают свойством симметрии. Таким образом, в общем случае динамическая жё- 
сткость процесса резания формирует силы, не обладающие потенциальными свойствами. Вместе с 
тем, существует режим (на рис.5 этот режим показан пунктирной линией), на котором матрицы 
динамической жёсткости характеризуются симметричными свойствами. В этом случае на вирту- 
альных перемещениях позиционные силы не совершают работу. Однако в этом случае матрица 
динамической жёсткости становится отрицательно определённой, то есть рассматриваемая точка 
равновесия является неустойчивой. 

При обработке резанием инструментами, обладающими значительной изгибной жёстко- 

стью, и во всех рассмотренных случаях матрицы динамической жёсткости принципиально зависят 
от геометрии инструмента. От геометрии зависит и потенциальность позиционных сил, форми- 
руемых процессом резания. В заключение отметим, что указанные свойства непотенциальности 
матриц динамической жёсткости процесса резания, а также перераспределения элементов матриц 
динамической жёсткости зависят и от скорости резания. 
Заключение. Обобщая приведённый материал, необходимо отметить, что параметры динамиче- 
ской характеристики процесса резания не являются инвариантными к геометрическим парамет- 
рам режущего инструмента и упругим свойствам подсистемы инструмента. Варьирование геомет- 
рии режущего инструмента и ориентации эллипсов жёсткости подсистемы режущего инструмента 
вызывает изменения матрицы жёсткости системы резания. В частности, на основании изменения 
только геометрии режущего инструмента можно существенно влиять на смещение точки равнове- 
сия системы резания, которое определяет упругие деформационные смещения инструмента отно- 
сительно заготовки, то есть на геометрию формируемой при резании детали. При этом вариации 
как геометрических параметров инструмента, так и эллипсов жёсткости подсистемы режущего 
инструмента изменяют условия самовозбуждения динамической системы резания. Это обусловле- 
но влиянием циркуляционных сил, формируемых процессом резания. Кроме этого, использование 
традиционно рассматриваемого коэффициента резания не позволяет раскрыть сложные преобра- 
зования силовых установившихся реакций со стороны процесса резания на деформационные 
свойства системы. 
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